Multiple Gaps and Superfluid Density from Interband Pairing in Iron Oxypnictides 
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We study a four-band model for the iron oxypnictides, in which the superconducting properties 
are assumed to be determined by the interband coupling between hole-like and electron-like Fermi 
sheets. We show that reasonable parameters can account for the angle-resolved photoemission spec- 
tra showing multiple gaps in Bai_ I K I Fe2As2, and for the temperature dependence of the superfluid 
density. At the same time, the zero-temperature value of the superfluid density shows a conventional 
scaling with the number of carriers. 
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The discovery of superconductivity in the family of 
compounds characterized by Fe-As layers^ is intriguing 
and challenging. Moreover, at the moment they are 
the only materials with T c larger than 40K beside the 
cuprates. While some properties of these materials are 
reminiscent of the cuprates (superconductivity induced 
by doping a magnetic parent compound, layered basic 
structure), the differences are also remarkable. 

Among the peculiarities of the iron oxypnictides not 
found in cuprates we mention the magnetic nature of Fe 
and the multi-band character of the band structure. The 
effect of correlations seems to be weaker than for the 
cuprates 2 -, even if also in the latter case the correlation 
strength may be smaller than usually believed 3 -. At the 
same time, the calculated electron-phonon coupling can 
not account for the high critical temperatures*, and some 
more exotic mechanism is expected to be realized. 

The initial efforts in the field focused on ROi-^F^FeAs 
and ROi-^FeAs (R being a rare earth atom), in which 
electrons are doped in the parent ROFeAs. Even more 
recently, a new direction opened thanks to the synthesis 
of Bai_ 2; K a; Fe2As2 and Sri_ a; K :c Fe2As2, doped by holes. 
The great advantage of this class of materials is the avail- 
ability of very large single crystals, which allows for a 
much more complete and reliable experimental analysis. 

Density functional theory (DFT) calculations of the 
band structure of oxypnictides show that the Fermi sur- 
face (FS) is formed by disconnected hole-like pockets cen- 
tered around the T point, and electron- like pockets cen- 
tered around the M points of the folded Brillouin zone 
of the FeAs planes. All the pockets essentially arise from 
the five d-orbitals of Fe, three pockets belonging to the 
hole-like manifold and two forming the electron- like one 5 . 
Such a FS topology, along with the magnetic character of 
the parent compound, motivated one of the first proposal, 
that assigned to the interband coupling between the two 
manifold of FS's the role of driving force behind the su- 
perconducting phenomenon 5 '^.. From a microscopic point 
of view, the origin of the interband pairing can be asso- 
ciated to spin fluctuations. Indeed, the magnetic order- 
ing vector of the parent compounds roughly connects the 



two sets of hole and electron FS's, leading to an (approx- 
imate) nesting, eventually giving rise to an effective pair- 
ing interaction. On the other hand, attractive phonon- 
induced intraband pairing (or the pairing between FS's 
of the same nature) is likely to be relatively small4 In 
this perspective, different nesting properties between the 
various FS sheets can be expected to lead to different 
strengths of the pairing interaction in the various bands. 

Recent angle-resolved photoemission (ARPES) results 
for undopedi and hole-doped x = 0.4 BaFe2As2^ clearly 
show that this can indeed be the case. From a broad per- 
spective, ARPES confirms the prediction of DFT about 
the existence of distinct electron and hole- like pockets. 
However, only two out of the three hole-like pockets may 
be resolved around the T point, with considerable dif- 
ferent areas, and a significant mass renormalization is 
found in all the bands with respect to DFT^ The su- 
perconducting (SC) gap opens on all the FS sheets at 
the same temperature, with a nodeless structure essen- 
tially compatible with simple s-wave symmetry. More 
remarkably, the gap at the smallest hole-like FS and the 
one at the electron-like sheets are almost identical, while 
the largest hole-like FS has half of this gap. As we shall 
see below, this effect is highly non trivial, and requires a 
proper analysis of the multiband model and of the role 
of interband interactions. The experimental situation is 
less clear for the electron-doped materials. A nodeless SC 
gap at the T-FS has been reported for NdFeAsOo.gF .i— , 
but no information is available to our knowledge about 
the gap in the electron-like pocket (s). 

Somehow related to the issue of the multiband struc- 
ture is the behavior of the superfluid density p s , both at 
zero and finite temperature. The T = BCS value of p s 
is indeed related to the filling and effective masses of the 
various bands, while its temperature dependence probes 
via quasiparticle excitations the presence of distinct SC 
gaps. The aim of this paper is to study a simple model 
which can capture the main effects of the multiple-gaps 
opening in BaFe2As2 superconductors. We will show that 
a reasonable description of ARPES results and superfluid 
density can be achieved within a four-band model with 
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realistic parameters, and that a doping-dependent renor- 
malization of the band structure should be envisaged. 
Finally, we show that the calculated T — value of the 
superfluid density is in reasonable agreement with the 
experiments ) 11 ' 12 ! 13 ' 14 ! 15 in striking contrast to cuprates, 
where the anomalous scaling of p s with doping can only 
be accounted for by including correlation effect o 16 ! 17 . 




FIG. 1: (Color Online) Schematic of the multiband model 
used in this work. The two hole bands a, /3 are centered 
around the T point, the electron bands around the M point. 

We consider a four-band model with parabolic disper- 
sion, £i(k) = e l Q ± ti(k — ki) 2 , with the plus (minus) sign 
for the electron (hole) bands and U hopping parameter 
(we will use, unless explicitly stated, units h = a = 1, a 
being the lattice spacing) . For the two hole bands around 
the r point, labeled as a and (3 (according to the nota- 
tion of Ref. [8]), kj =0, while for the two electron bands 
7 around the M points k^ = (tt,tt). As we mentioned 
above, the third hole pocket predicted by DFT is not 
seen around the T point, probably because highly degen- 
erate with the a one. Moreover, the two hole pockets are 
slightly split in the undoped systemj 7 - Thus, we choose 
the reference energy such that — 0, while Eq and eg 
are in general different from zero, see Fig. [T] For sim- 
plicity, we assume that the two electron bands are equal 
and isotropic around the M points, even though both 
DFT and experiments show that they are slightly elon- 
gated along the k x and k y directions. Within a BCS-like 
approach, the pairing Hamiltonian can be written as: 

H = e *(k)c+ k(T c l!k(7 + ^,i c tkT c ?-ki c J^k'iCi,k'T,(l) 

i,k<r 

k,k' 

where cf k creates an electron in the i-th band (i = 
a, (3, 7), and as we discussed above we neglect intraband 
pairing. Moreover, the similar size of the FS for the a 
and 7 bands suggests that the largest pairing acts be- 
tween these two FS sheets, which are strongly nested, 
while the f3 and 7 bands are less coupled. These condi- 
tions are implemented in by taking V a p — 0, Vp~ ( — A 
and V ai = A > A. As it has been noticed in Ref.-£, for 
inter-band mechanisms it is not necessary that the inter- 
action is attractive (at least at low energy) , as it happens 
for intra-band mechanisms. However, if A, A > one sees 
that the order parameters must have different sign on dif- 
ferent kind of FS sheets, i.e. signA Qi( g = — signA 7 . The 
resulting set of self-consistent equations for the absolute 



values of the gaps is a straightforward extension of the 
usual BCS equations: 

A a = AA 7 2II 7 (2) 
A/3 = AA 7 2II 7 (3) 
A 7 = AA Q Il a + AA^IIg (4) 

where A^ is the gap in the z-th band, whose disper- 
sion is now = y / |f(k) + A 2 , with & = e< - fi. 
The IF are the particle-particle bubbles, given by IF = 
N l J^°d^(l/E l )ta,nh(E l /2T) 1 where the pairing is only 
effective for states in a range |£| < ojq around the Fermi 
level and N, — 1/Airbi is the density of states (DOS) of 
each band. Here we use ujo — 20 meV, i.e. of the or- 
der of the smallest Fermi energy in the different sheets. 
Along with Eqs. ©-([1]) we must solve self-consistently 
the equation for the particle number n = n a + np + 2n 7 , 
where n = 4 is the half filling. At the present stage the 
agreement between the experimental dispersion and the 
DFT bandstructure is hardly satisfactory^, therefore we 
extract the band parameters from the experimental data 
of Ref. 7 for undoped BaFe2As2. Here one has t a = 90 
meV, tp — 39 meV and t 7 = 45 meV, with a band split- 
ting £q = —10 meV and Eq = —45 meV. However, in a 
rigid-band picture, with these parameters at the doping 
x = 0.4 the electron pockets would be empty. Moreover, 
as we shall see below in Fig. ^jp, with these parameters 
the gap values in the a and 7 bands would be quite dif- 
ferent, in contrast to what measured in Ref.~. We thus 
proceed by assuming a band renormalization with respect 
to half-filling to obtain A Q ~ A 7 , and A a = 2Ap from 
Eqs. ©-(HI), where the gap values depend only on the 
DOS and on the coupling. The latter condition is always 
realized for A = A/ 2. The former condition instead can 
be obtained with slightly renormalized bands, i.e. we use 
t a = 65 meV, tp = 25 meV and t 1 = 60 meV. Notice 
that if a third hole pocket of the same size of the a one is 
present, even if not explicitly observed by ARPES, this 
would just require to slightly change the parameters to 
compensate the substitution of n Q — > 2n Q in Eq. {H . 

Once the DOS are fixed by the gap values, we tune the 
band splitting according to the FS areas Ai measured 
in RefA In particular, it is observed that A a /Ap = 
2/9. Notice that if a and /3 bands were degenerate (or 
nearly degenerate) at the T point this A a /Ap ratio would 
require a very large anisotropy between the hole bands, 
tp = i a (2/9), which in turn would not allow the condition 
A Q ~ A 7 to be fulfilled. Thus, to justify the ratios of 
the FS areas we should assume a splitting of the hole 
bands at the T point, even larger than what measured 
at half-filling 7 -. This enhancement of the splitting can be 
associated to a different Hartree shift of the two bands, 
arising from a particle-hole decoupling of the interband 
interaction. Imposing the FS ratios suggested in Ref4 
we determine Cq = —71 meV and Eq = —25 meV (so that 
ji ~ — 60 meV). In Fig. [2^ we show the results for the 
gaps using A = 260 meV, which at low T are in good 
agreement with the experimental results of Ref£, even if 
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FIG. 2: (Color Online) (a): Temperature dependence of the 
SC gaps in the hole and electron bands. The parameter val- 
ues are optimized to obtain the same gap value in the a and 
7 bands (values in eV). Moreover, we used A = A/2 to have 
A^ = A a /2. (b): SC gaps for the hopping parameters esti- 
mated from ARPES measurements of Ref.— at half-filling. 

a larger T c is found here because fluctuation effects, that 
are expected to be important in 2D, are not included. We 
point out that in a multiband system dominated by off- 
diagonal (interband) coupling, deviations from the BCS 
value of 2A/T C in each band can be obtained without 
invoking exotic mechanisms. In Fig. [5^ we find 2A^/T C = 
2 and 2A a /T c = 4, and further enhancement could be 
expected after the inclusion of fluctuation effects. 

When multiple gaps are present, one may expect a sig- 
nature on the temperature behavior of the superfluid den- 
sity p s . For a generic lattice dispersion ek, the superfluid 
density at BCS level is given by: 

[<Pk fde k \ 2 df(E w ) 

+ 2 J ~^\dk-J ~mT' (5) 

where fix) = [1 + exp(x/T)] _1 is the Fermi function. 
Let us first consider the electron bands 7: since they 
are almost empty, we can approximate the full disper- 
sion ek with the parabolic band e 7 introduced above, 
so that the previous equation reduces to p], = 2i 7 n 7 + 
4i 7 7V 7 J deef'iEj). As T — > the quasiparticle excita- 
tions vanish, due to the derivative f'(Ej), and the usual 
result pi — 2i 7 n 7 is recovered. For hole bands, a simi- 
lar expression can be used after a particle-hole transfor- 
mation, which insures that only the hole states at the 
top of the band contribute to the superfluid density, i.e. 
p a s = 2t a (2-n a ) at T = 0. 

The results for the temperature dependence of the su- 
perfluid density for the same parameters of Fig. [2^ are 
shown in Fig. [3J where we plot the total p s together with 
the contributions from the single bands. In the (3 band, 
which has the smallest gap, the clear signature of the 
two different order parameters can be seen most clearly, 
while the two other bands have a more standard behav- 
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FIG. 3: (Color Online) Temperature dependence of the su- 
perfluid density in units of Eq. (|6]) for the various bands. 



ior. As we observed above, for the /3-band 2Ap/T c = 2, 
i.e. it is smaller than the BCS value. This means that 
A(T) is finite and sizable in a regime of temperature 
where it would be zero if not coupled to another band. 
As a consequence, at low T p@ decreases quite rapidly, 
as it would vanish at a temperature proportional to A^, 
while at higher temperature the depletion softens in a tail 
which closes at the higher T c determined by the inter- 
band coupling. The possibility to observe a clear change 
of curvature in the total superfluid response p' ot depends 
on the relative filling of the band with the smaller gap. 
In Fig. [3J we show that the presence of the multiple gaps 
leads to deviations of p s ot (T) with respect to the single- 
band case, but with no upward curvature as in p@(T). 
We notice that pi"* strongly resembles the one reported 
recently in Refji 1 - in a Bao.55Ko.4sFe2As2 compound with 
similar doping of the one studied by ARPES in Ref. 8 . 

More intriguing is instead the temperature dependence 
observed in several electron-doped LaOi-^FaFeA s 12 ^ 3 
and NdOi-rcFzFeAs^ compounds, where the measured 
p s resembles the one we find for the (3 band. This effect 
is particularly pronounced in LaOi-^F^FeAs for x > 0.1, 
and suggests that the smaller gap is formed on the elec- 
tron part of the FS, which is expected to give the larger 
contribution to the superfluid response in electron-doped 
compounds. Even though some indications of multiple 
gaps in LaOi-^F^FeA*^ already exist, a clear experi- 
mental observation of the gap opening on the different 
FS sheets is still lacking. 

We finally comment on the T = value of the su- 
perfluid density, which is proportional to the Fermi en- 
ergy within our BCS scheme, and it is therefore of 
the order of several hundreds of K. This energy scale 
is comparable to the values currently available from 
penetration-depth measurements in both hole-dope d 11 ! 15 
and electron-dope d 12 ' 13 i 14 oxypnictides. To make a di- 
rect comparison with the experiments, we recall that in 
two dimensions the quantity J s = fr?n 2d jm is an energy, 
where n 2 s d is a two-dimensional density of superfluid car- 
riers. By considering each plane as the basic 2D unit, 
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one can relate n 2 and J s to the measured penetration 
depth A as: 

_ h 2 nl d _ h 2 c 2 d _ d[A] 
Js ~ ~ 47re 2 A 2 " ~ 2M \ 2 {{p,m) 2 ] K (6) 

where d is the interlayer spacing, which ranges between 
4.3-6.6 A in oxypnictides. The penetration depth of a 
Bao.6Ko.4Fe 2 As 2 sample has been estimated in Ref. 15 
from infrared measurements as A = 0.208 pm. This gives 
according to Eq. |6|) J s — 380 K, which is not much lower 
that the value we estimated above assuming the clean 
limit. Moreover, since the optical-conductivity measure- 
ments of Reffi^ suggest that the sample is in the dirty 
limit, we would expect that only a fraction of the p s esti- 
mated from the (renormalized) band structure will actu- 
ally condense in the SC phase. We can then conclude that 
the measured superfluid-density values are quite consis- 
tent with band-structure calculations. For this reason, 
the observed similarity of the relation of 1/A 2 vs T c in 
oxypnictides and cuprates superconductors^ 2 - should not 
be taken as a signature of similar unconventional behav- 
ior in the two classes of materials. Indeed, in cuprate 
superconductors the number of carriers is large, because 
one is doping the system with respect to half-filling, 
where the superfluid density is expected to attain its 
maximum value in a hypothetical uncorrelated system. 
In other words, one would expect p s ~ (1 — x)t, while 



one measures values of order p s ~ xt, where x is the 
doping and t the typical hopping parameter. For this 
reason, correlation effects as the proximity to a Mott 
insulator— 6 or enhanced phase fluctuations^ 7 , have been 
invoked as possible explanation of the low superfluid den- 
sity in cuprates. Instead, in oxypnictides the system 
has a low superfluid density because an almost full hole 
band has a vanishing contribution to p s , and one expects 
p s ~ xt which is indeed near to the experimental values. 

In conclusion, we have presented calculations in a four- 
band model built on the basis of the experimental results 
on BaFe2As2, assuming that the pairing arises from the 
interaction between hole-like and electron-like Fermi sur- 
faces, with different strengths due to different nesting 
properties of the various FS sheets. Within our multi- 
band BCS approach we determine parameters that ac- 
count for the multiple gaps recently observed by Ding 
et al. s . Using the same parameters, we compute the su- 
perfluid density, which does not show striking signatures 
of the multiple gaps, in agreement with experiments. 
The same agreement holds for its zero-temperature value, 
which scales with the small density of carriers. This has 
to be contrasted to the case of cuprates, where the small- 
ness of p s is associated to correlation effects. 

We acknowledge R. Arita, L. Boeri, X. Dai, and 
M. Grilli for useful discussions. 
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